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Abstract 

We implement an algorithm which is aimed to reduce the number 
of basis states spanning the Hilbert space of quantum many-body 
systems. We test the procedure by working out and analyzing the 
spectral properties of strongly correlated and frustrated quantum spin 
systems. The role and importance of symmetries are investigated. 
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1 Introduction. 

Most microscopic many-body quantum systems are subject to strong inter- 
actions which act between their constituents. The description of some sys- 
tems can be tackled by means of perturbation theory. This may be the case 
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when it is possible to introduce a mean-field concept which is able to include 
quantitatively the main part of the interactions, leaving a remaining residual 
contribution which acts as a more or less small perturbation. 

Often such an approach does not lead to sensible results, in particular 
when one is dealing with realistic quantum spin systems. Non-perturbative 
techniques are needed. During the last decades a considerable amount of 
procedures relying on the renormalization group concept introduced by Wil- 
son pQ have been proposed and tested. Some of them are specifically de- 
vised for quantum spin systems, like the Real Space Renormalization Group 
(RSRG) H E] and the Density Matrix Renormalization Group (DMRG) 

In many cases the study of the spectral properties of quantum systems is 
obtained through the diagonalization of a many-body Hamiltonian in Hilbert 
space spanned by a complete, in general infinite or at least very large set of ba- 
sis states although the information of interest is restricted to the knowledge 
of a few low energy states generally characterized by collective properties. 
Consequently it is necessary to manipulate very large matrices in order to 
extract a reduced quantity of informations. 

Recently we proposed a non-perturbative approach which tackles this 
question [2J. The procedure consists of an algorithm which implements a 
step by step reduction of the size of Hilbert space by means of a projection 
technique. It relies on the renormalization concept following in spirit for- 
mer work based on this concept [91 [10l [11] . Since the reduction procedure 
does not act in ordinary or momentum space but in the zero- dimensional 
Hilbert space like in the procedure developed in ref. [TS], it is in principle 
applicable to all types of microscopic quantum systems, in contradistinction 
with methods like the DMRG procedure which works in ordinary space and 
is specifically applicable to Id quantum lattice systems. 

In the present work we implement this algorithm as a preliminary test of 
the practical efficiency and accuracy of the method when applied to strongly 
interacting systems which cannot be treated by means of perturbative meth- 
ods. Second we want to see how far it is able to deliver physical information 
about the properties of the many-body systems it is aimed to describe. We 
use quantum spin systems as first test probes. 
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The outline of the paper is the following. In section 2 we recall the 
essential steps leading to the derivation of the equation which governs the 
evolution of the coupling strength attached to the interaction. The reduction 
formalism is universal in the sense that it works for any kind of many-body 
quantum system. Section 3 is devoted to the application of the algorithm 
to frustrated quantum spin ladders with two legs and one spin per site. We 
analyze the outcome of the applied algorithm on systems of different sizes 
and characterized by different coupling strengths by means of numerical ex- 
amples, with bases of states developed in the SU(2) and S'0(4)-symmetry 
scheme. General conclusions and further planned investigations and devel- 
opments are drawn in section 4. 

2 The reduction algorithm. 

2.1 General concept: the space reduction procedure. 

We consider a system of quantum objects (particles, spins) which are char- 
acterized by a discrete spectrum. The system is governed by a Hamiltonian 



|<7i , #2 j ' ' ' > 9p i— *■ an d acts in a Hilbert space of dimension 

N. The spectrum is obtained from 



I (g( N ^)), i = 1, • • • ,Nj depend on the set of coupling constants 

{g( N ^}- If the relevant quantities of interest are for instance M eigenval- 
ues out of the original set it makes sense to try to define a new effective 
Hamiltonian H^ M \g^) whose eigenvalues reproduce the M selected states 
and verifies 





(1) 



(2) 



with the constraints 



W M) ) = W N) ) 



(3) 
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for % = 1,...,M. If this can be realized Eq. (3) implies a relation between 
the coupling constants in the original and reduced space 

n (M) _ , f (2V) (N) (AT)n (A s 

with k = l,...,p. The effective Hamiltonian H^ M \g^ M ^) may not be rigor- 
ously derivable from H^ N \ It should be constructed so that it optimizes the 
overlap between the original and reduced set of eigenstates. We show next 
how this space reduction may be implemented in practice. 



2.2 Reduction procedure and renormalization of the 
coupling strengths. 

We sketch the procedure which leads from Eq. flTJ to Eq. (T5]). Details can be 
found in ref. [2J. 

We consider a system described by a Hamiltonian depending on a unique 
coupling strength g which can be written as a sum of two terms 

H = H + gH 1 (5) 

The Hilbert space HS N > of dimension N is spanned by an a priori arbitrary 
set of basis states % — 1, • • • , iV}. An eigenvector ) can be written 

as 

N 

\< N) ) = J2al 1 I\9 {N) m) (6) 
i=i 

where the amplitudes {a^ig^)} depend on the value of g in 

Using the Feshbach formalism [12J the Hilbert space may be decomposed 
into subspaces by means of the projection operators P and Q, 

H (N) = pH (N) + QH (N) (7) 

In practice the subspace P7i^ is chosen to be of dimension dim PJi^' = 
N — 1 by elimination of one basis state. The projected eigenvector P\^/[ N ^} 
obeys the Schroedinger equation 

H eff (x[ N) )P\^) = . (8) 
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where H e ff(\\ ) is the effective Hamiltonian which operates in the sub- 
space PT~L^ N \ It depends on the eigenvalue Aj which is the eigenenergy 
corresponding to N ') in the initial space TiS N \ The coupling which 
characterizes the Hamiltonian H^ N > in TC^ N > is now aimed to be changed into 
such a way that the eigenvalue in the new space TiS N ^ is the same 
as the one in the complete space 

\[ N - 1] = \[ N) (9) 

The determination of g^ N ~^ by means of the constraint expressed by Eq. (Q 
is the central point of the procedure. In practice the reduction of the vector 
space from N to N — 1 results in a renormalization of the coupling constant 
from to g^ N ~^ preserving the physical eigenenergy \[ N \ 

In the sequel P^^) is chosen to be the ground state eigenvector and 
Aj^ = \[ N ^ = Ai the corresponding eigenenergy. In ref. |2j it is shown how 
gC/v-i) can De obtained as a solution of an algebraic equation of the second 
degree, see Eq. (17). The reduction procedure is iterated in a step by step 
decrease of the dimensions of the vector space, N\—>-N— 1*-^N — 2*-^... 
leading at each step Ho a coupling strength g( N ~ k ) which can be given as 
the solution of a flow equation in a continuum limit description of the Hilbert 
space [2j [2TJ . The procedure can be generalized to Hamiltonians depending 
on several coupling constants. 



2.3 Reduction algorithm. 

1— Consider a quantum system described by an Hamiltonian which 
acts in an N- dimensional Hilbert space. 

2— Compute the matrix elements of the Hamiltonian matrix in a def- 
inite basis of states = 1, . . . , N}. The diagonal matrix elements 
{cj = (<&i\H N \<$>i)} are arranged in increasing order with respect to the {ej}. 

3— Use the Lanczos technique to determine A^ and |^i (<7^))- A^ may 
be chosen as the experimental value Ai if known. 

4— Fix g^ N ~^ as described in section 2.2. Take the solution of the algebraic 
second order equation closest to g( N ) (Eq. (17) in [2]). 

5— Construct if^" 1 ) = Hq-\- g^ N ~ 1 ' Hi by elimination of the matrix elements 
of H^ N > involving the state |3>jv). 
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6— Repeat the procedures 2, 3, 4 and 5 by fixing at each step k X\ = 
X[ N) = A, 

7— The iterations may be stopped at iV = N min corresponding to the limit 
of space dimensions for which the spectrum gets unstable. 

2.4 Preliminary remarks. 

• The procedure is aimed to generate the energies and other physical 
properties of the ground state and low-energy excited states of strongly 
interacting systems. 

• The implementation of the reduction procedure asks for the knowledge 
of Ai and the corresponding eigenvector at each step k of the 
reduction process. The eigenvalue Ai is chosen as the physical ground 
state energy of the system. Eigenvalue and eigenvector can be obtained 
by means of the Lanczos algorithm [T5l [T6l [8] which is particularly well 
adapted to very large vector space dimensions. This algorithm is used 
here in order to fix A i and ). 

• The process does not guarantee a rigorous stability of the eigenvalue Ai. 
Indeed one notices that |\|/ 1 Ar_fc ~ 1 * ) ) which is the eigenvector in the space 
H (N-k-i) and the p ro j ectec i state P\^/[ N ~ k) ) of into U^'^ 
may differ from each other. As a consequence it may not be possible 
to keep A^ 1 ^ rigorously equal to A x = \\. In practice the degree 
of accuracy depends on the relative size of the eliminated amplitudes 
{dyjig^)}- This point will be tested by means of numerical estima- 
tions and further discussed below. 

• The algorithm and different points which have been quoted above will 
now be developed in applications of the procedure to explicit models, 
here frustrated spin ladders. The applicability of the algorithm will be 
tested in different symmetry schemes. 

• The Hamiltonians of the considered ladder systems are characterized by 
a fixed total magnetic magnetization M tot . We shall work in subspaces 
which correspond to fixed M tot . The total spin S to t is also a good 
quantum number which defines smaller subspaces for fixed M tot . We 
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1 2 3 4 5 / L 

Figure 1: Top: the original spin ladder. The coupling strengths are indicated 
as given in the text. Bottom: The ladder in the SO(4) representation. See 
the text. 



shall not introduce them here because projection procedures on Stot are 
time consuming. Furthermore we want to test the algorithm in large 
enough spaces although not necessarily the largest possible ones in this 
preliminary tests considered here. 



3 Application to frustrated two-leg quantum 
spin ladders. 

3.1 The model 

3.1.1 SU(2)-symmetry framework. 

Consider spin- 1/2 ladders [12], UJ] described by Hamiltonians of the following 
type and shown in Fig. 1. 

L 

H (s, s ) = j t y^ 8ii8ia + j l s il s jl + j t s i2 s h + Jicy^s^sj,, (10) 

i=l <ij> <ij> (ij) 

{ij) 

The indices 1 or 2 label the spin 1/2 vector operators Si k acting on the 
sites i on both ends of a rung, in the second and third term i and j label 
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nearest neighbours, here j = i + 1 along the legs of the ladder. The fourth 
and fifth term correspond to diagonal interactions between sites located on 
different legs, j — i + 1. L is the number of sites on a leg (Fig. 1) where 
J\c = J2c = J c - The coupling strengths J t , Ji, J c are positive. 

As stated above the renormalization is restricted to a unique coupling 
strength, see Eq. It is implemented here by putting H = and = 
g( N *>Hi where = Jt and 

L 

Hi = ^ S h S i2 + Itl 5^ ( S h S jl + S i2 S 32) + 7c 5^ ( S h S j2 + S l2 S h) ■ (11) 
i=l <ij> <ij> 

where j t l — Ji/Jt, 7c = Jc/ Jt- These quantities are kept constant and 
g( N ) — j t w iU be subject to renormalizationin the reduction process. 

The basis of states {\§k) k = 1, . . . , N} is chosen as 

2L 

|$ fc ) = |l/2 mi, ...,1/2 m i; ...,l/2 m 2L , ^m, = M to4 = 0) 

i=l 

with {mi = +1/2,-1/2}. 

3.1.2 SO(4)-symmetry framework. 

The basis of states may be written in an SO(4)-symmetry scheme. Differ- 
ent choices of bases may induce a more or less efficient reduction procedure 
depending on the strength of the coupling constants J t , Ji, J c - This point is 
investigated here. 

By means of a spin rotation [TT1 [TS] 

Si^^Si + Ri) . (12) 

Si^^Si-Ri) . (13) 
the Hamiltonian Eq.(10) can be expressed in the form 
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H (s, R) = I - Ri) + Jx £ S t S 3 + J 2 £ R*R> (14) 

i=l <ij> <ij> 

The structure of the corresponding system is shown in the lower part of Fig. 
1. Here J\ = (Jj + J c )/2, J 2 = (J; — J c )/2 and as before J\ c = J 2c = J c - The 
components S^ + \ S^~\ and R\ + \ Ri~\ R^ of the vector operators Si and 
Ri are the 5*0(4) group generators and < ij > denotes nearest neighbour 
indices 

Sf > = V2(XJ U)(10) +^ 10)(1_1) ) = 



s-w = x (11)(11) - x (1_1)(1_1) 



i^ 2) = _(x (10){00) + x (00){10) ) 

where 

X (WW) = | 5 . M . )(W | 
In this framework the states {\SiMi)} are defined as 

iS'iMi) = ^ (1/2 mi 1/2 771213^)11/2 m^ll/2 m 2 )i 

m 1 ,m,2 

along a rung are coupled to Si — or Si — 1. Spectra are constructed in this 
representation as well as in the SU (2) representation and the basis of states 
{|$fc)} takes the form 

L 

|$ fc ) = l^iMi, SiMi, S L M L , Mi = M tot = 0) 

i=l 



9 



3.2 Test observables 

In order to quantify the accuracy of the procedure we introduce different 
test quantities in order to estimate quantitatively deviations between ground 
state and low excited state energies in Hilbert spaces of different dimensions. 
The stability of low-lying states can be estimated by means of 

/ (AO _ (N-k)\ 

p(i) = \¥i J| L\ x 100 with i=l,..., 4 (15) 

where e\ N ~ k ^ = \\ N ~ k ^ /2L corresponds to the energy per site at the ith phys- 
ical state starting from the ground state at the fcth iteration in Hilbert space. 
These quantities provide a percentage of loss of accuracy of the eigenenergies 
in the different reduced spaces. 

A global characterization of the ground state wavefunction in different 
symmetry schemes can also be given by the entropy per site in a space of 
dimension n 

n 

s = -^j2 p ^ with ^ = i<^ B) i*l B) >i 2 = i«& ) r ( 16 ) 

1=1 

which works as a global measure of the distribution of the amplitudes {a^} 
in the physical ground state. 



3.3 Spectra in the SU(2)-symmetry framework. 

We apply the reduction algorithm to ladders with two legs, different num- 
bers of sites and different values of the coupling strengths. Results obtained 
with an SU (2)-symmetry basis of states are shown in Figs. (2 - 4) and Fig. 6. 

3.3.1 First case: L= 6, J t =15, J t =5, J c =3 

We choose the basis states in the framework of the M-scheme corresponding 
to subspaces with fixed values of the projection of the spin of the {|<E»j)}, 
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M tot = 0. 



In the present case J t >Ji, J c - The dimension of the subspace is reduced 
step by step as explained above starting from iV = 924. As stated in section 
2.3 the basis states {|$?,)} are ordered with increasing energy of their diago- 
nal matrix elements {q} and eliminated starting from the state with largest 
energy, e N . 

As seen in Fig. (j2k) the ground state of the system stays stable down to 
n ~ 50 where n is the dimension of the reduced space. The coupling constant 
J t does not move either up to n ~ 300. Figs.([2h-b) show the evolution of the 
first excited states which follows the same trend as the ground state. Devia- 
tions from their initial value at iV = 924 can be seen in Fig.(T5b-d) where the 
p(i)'s defined above represent these deviations in terms of percentages. 

For n < 50 the spectrum gets unstable, the renormalization of the cou- 
pling constant can no longer correct for the energy of the lowest state. Indeed 
the coupling constant J t increases drastically as seen in Fig.([2fe). The reason 
for this behaviour can be found in the fact that at this stage the algorithm 
eliminates states which have an essential component in the state of lowest 
energy. The same message can be read on Fig.([2]F), the drop in the entropy 
per site s is due to the elimination of sizable amplitudes {an}. 

3.3.2 Second case: L= 6, J t —5.5, J;=5, J c =3 

Contrary to the former case the coupling constant Jt along rungs is now of 
the same strength as Ji, J c . Results are shown in Fig. ([3]). The lowest energy 
state is now stable down to n ~ 100. This is also reflected in the behaviour 
of the excited states which move appreciably for n < 200. Fig.([3fe) shows 
that the coupling constant J t starts to increase sharply between n = 300 and 
n = 200. It is able to stabilize the excited states up to about n = 200 and 
the ground state up to n = 70. The instability for n < 70 reflects in the 
evolution of the p(«)'s, Figs.([3b-d) which get of the order of a few percents. 
The entropy Fig.([3]F) follows the same trend. 

Comparing the two cases above and particularly the entropies Fig.([3]F) 
and Fig.(T5]F) one sees that the stronger J t the more the amplitude strength of 
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the ground state wavefunction is concentrated in a smaller number of basis 
state components. The elimination of sizable components of the wavefunction 
leads to deviations which can be controlled up to a certain limit by means 
of the renormalization of Jt- One sees that large values of Jt favour a low 
number of significative components in the low energy part of the spectrum 
in a SU(2) symmetry framework. 

A confirmation of this trend can be observed in Figs.(jlK-f) where Jt = 2.5. 
The rates of destabilisation of the excited states are higher than in the for- 
mer cases as it can be seen in Figs.(JIb-d). This point is also reflected in the 
behaviour of the entropy s which is larger than in the former case for n = N 
and decreases more rapidly with decreasing n, Fig.(J3f). 

3.3.3 Third case: L= 8, J<=15, J,=5, J c =3 

For L = 8 the Hilbert space is spanned by iV = 12870 basis states with 
M to t = 0. The results are shown in Figs.([6K-f). The stability of the spectrum 
with decreasing space dimension is stronger than the stability observed for 
L = 6. Indeed if n/N defines the ratio of the number of states in the reduced 
space over the total number of states one finds p{l) ~ 0.8% and p{2) ~ 0.8% 
when n/N ~ 0.07 for L = 6. For L = 8 p(l) ~ 0.8% when n/N ~ 0.007 and 
p(2) ~ 0.5% when n/N ~ 0.02. This shows a sizable improvement in the 
stability of the spectrum, at least in the specific domain where the coupling 
strength J t is large compared to the others. 

The evolution of the spectrum and its stability with decreasing J t follows 
the same trend as in the case where L = 6. 



3.3.4 Remarks 

In Fig.(TSk) it is seen that the ground state shows "bunches" of energy fluc- 
tuations. The peaks are intermittent, they appear and disappear during the 
space dimension reduction process. They are small in the case where Jt = 15 
but can grow with decreasing J t as it can be observed for J t = 2.5. The 
subsequent stabilization of the ground state energy following such a bunch 
shows the effectiveness of the coupling constant renormalization which acts 
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in a progressively reduced and hence incomplete basis of states. 

These bunches of fluctuations are correlated with the change of the num- 
ber of relevant amplitudes (i.e. amplitudes larger than some value e as ex- 
plained in the caption of Fig.©) during the reduction process. 

Consider first the case where J t > Ji, J c - One notices in the caption of 
Fig.©i) that up to n ~ 300 the number of relevant amplitudes defined in 
Fig.© stays stable like the ratios {p(i)} in Figs.([2b-d). For 158 < n < 300 
these ratios change quickly. A bunch of fluctuations appears in this domain 
of values of n as seen in Figs. ([2b-d) and correspondingly the number of rel- 
evant amplitudes decreases steeply. For 60 < n < 158 the ratios {p{i)} stay 
again stable as well as the number of relevant amplitudes. The {p(i)} in 
Fig.(T5b-d) almost decrease back to their initial values. The same explanation 
is valid for L = 8 in Fig.©. The analysis shows that these bunches of fluctu- 
ations signal the local elimination of relevant contributions of basis states to 
the physical states in the spectrum. The stabilization of the spectra which 
follows during the elimination process shows that renormalization is able to 
cure these effects. 

In the case where J t < Ji,J c shown in Fig.©)) the relevant and irrel- 
evant amplitudes move continuously during the reduction process and the 
corresponding {p{i)} do no longer decrease to the values they showed before 
the appearance of the bunch of energy fluctuations as seen in Figs.©-d). It 
signals the fact that the coupling renormalization is no longer able to com- 
pensate for the reduction of the Hilbert space dimensions. 



3.4 Spectra in the SO(4)-symmetry framework 

The reduction algorithm is now applied to the system described by the Hamil- 
tonian 

H (S,R) 

given by Eq. (I14p with a basis of states written in the SO (A) 
symmetry framework. Like above we consider two cases corresponding to 
large and small values of Jt relative to the strengths of the other coupling 
parameters. 
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3.4.1 Reduction test for L = 6, J t = 15, 2.5 and J;=5, J c =3 

Figs.© show the behaviour of the spectrum for a system of size L = 6. A 
large value of Jt, (Jt = 15), favours the dimer structure along rungs in the 
lowest energy state and stabilizes the spectrum down to small Hilbert space 
dimensions. This effect is clearly seen in Fig.flTti), the ground state is very 
stable. The excited states are more affected, see Figs. (J7b), although they 
do not move significantly, Figs.([7b-d). The renormalization of the coupling 
strength J t starts to work for n ~ 50. 

The situation changes progressively with decreasing values of J t . Figs.([HD 
show the case where J t = 2.5. The ground state energy experiences sizable 
bunches of fluctuations like in the SU(2) scheme, but much stronger than 
in this last case. The same is true for the excited states which is reflected 
through all the quantities shown in Figs.©, in particular J t , Fig. ([Hb). The 
arguments used in the SU(2)-scheme about relevant and irrelevant ampli- 
tudes are also valid here. 

The result shows that the renormalization procedure is quite sensitive to 
the symmetry scheme chosen in Hilbert space. It is expected that essential 
components of the ground state wavefunction get eliminated early during 
the process when the rung coupling gets of the order of magnitude or smaller 
than the other coupling strengths. 



3.5 Summary 

The present results lead to two correlated remarks. The efficiency of the 
algorithm is different in different sectors of the coupling parameter space. 
In the case of the frustrated ladders considered here the algorithm is the 
more efficient the stronger the coupling between rung sites Jt- Second, this 
behaviour is strongly related to the symmetry representation in which the 
basis of states is defined. The SU(2) representation leads to a structure of 
the wavefunctions (i.e. the size of the amplitudes of the basis states) which is 
very different from the one obtained in the 5*0(4) representation. For large 
values of Jt the spectrum is more stable in the SO(4) scheme. For small 
values of J t the stability is better realized in the SU(2) scheme. Finally, in 
the regime where J t > Ji, J c , one observes that the reduction procedure is the 
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more efficient the closer J\ to J c . This effect can be understood and related 
to previous analytical work in the SO (4) framework |20j. 



4 Conclusions and outlook. 

In the present work we tested and analysed the outcome of an algorithm 
which aims to reduce the dimensions of the Hilbert space of states describing 
strongly interacting systems. The reduction is compensated by the renor- 
malization of the coupling strengths which enter the Hamiltonians of the 
systems. By construction the algorithm works in any space dimension and 
may be applied to the study of any microscopic iV-body quantum system. 
The robustness of the algorithm has been applied to frustrated quantum spin 
ladders. 

The analysis of the numerical results obtained in applications to quantum 
spin ladders leads to the following conclusions. 

• The stability of the low-lying states of the spectrum in the course of 
the reduction procedure depends on the relative values of the coupling 
strengths. The ladder favours a dimer structure along the rungs, i.e. 
stability is the better the larger the transverse coupling strength J t . 

• The evolution of the spectrum depends on the initial size of Hilbert 
space. The larger the initial space the larger the ratio between the 
initial number of states and the number of states corresponding to the 
limit of stability of the spectrum. 

• The efficiency of the reduction procedure depends on the symmetry 
frame in which the basis of states is defined. It appears clearly that the 
evolution of the spectrum described in an SU (2) scheme is significantly 
different from the evolution in an 5*0(4) scheme. This is again under- 
standable since different symmetry schemes partition Hilbert space in 
different ways and favour one or the other symmetry depending on the 
relative strengths of the coupling constants. 
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Local spectral instabilities appearing in the course of the reduction pro- 
cedure are correlated with the elimination of basis states with sizable 
amplitudes in the ground state wavefunction. One or another repre- 
sentation can be more efficient for a given set of coupling parameters 
because it leads to physical states in which the weight on the basis 
states is concentrated in a different number of components. This point 
is strongly related to the correlation between quantum entanglement 
and symmetry properties which are presently under intensive scrutiny, 
see f.i. [22] and refs. therein. 

Further points are worthwhile to be investigated: 

• In the present approach the sequential reduction of space dimensions 
followed an energy criterion. It might be judicious to classify the se- 
quence of states to be eliminated starting with those which have the 
smallest amplitude in the ground state wavefunction. The two proce- 
dures should however be correlated if not equivalent. 

• We expect to extend the study to systems of higher space dimensions, 
f.i. 2d. 

• The present approach relies on an algorithm which is able to recognize 
the existence of first and higher order critical points [2J. It is of in- 
terest to apply the algorithm in the neighbourhood of such points. Its 
behaviour could help to identify them. Work on this point is under way. 

• The algorithm can be extended to systems at finite temperature [21] 
and more than one coupling constant renormalization. 

The authors would like to thank Dr. A. Honecker for interesting discus- 
sions. One of us (T.K.) acknowledges the help of Drs. E. Caurier and F. 
Nowacki on technical aspects concerning the implementation of the Lanczos 
algorithm. 
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Figure 2: SU(2) — scheme. N is Hilbert space dimension. The {ei,i 
1,2,3,4} are the energies of the ground and excited states per site. L = 
sites along a leg. J t = 15, J\ = 5, J c = 3 
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Figure 3: SU(2) — scheme. N is Hilbert space dimension. The {e^, i 
1,2,3,4} are the energies of the ground and excited states per site. L = 
sites along a leg. J t = 5.5, Ji = 5, J c = 3 
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Figure 4: SU(2) — scheme. N is Hilbert space dimension. The {e^, i = 
1, 2, 3, 4} are the energies per site. L = 6 sites along a leg. J t = 2.5, Ji = 5, 

Jr. = 3 
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Figure 6: SU(2) — scheme. N is Hilbert space dimension. The {e^, i = 
1, 2, 3, 4} are the energies per site. L = 8 sites along a leg. J t = 15, J/ = 5, 
J, = 3 
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Figure 7: SO (4) — scheme. N is Hilbert space dimension. The {e^ 
1,2,3,4} are the energies per site. L = 6 sites along the chain. J t = 
Ji = 5, J c = 3 
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Abstract 

We implement an algorithm which is aimed to reduce the number 
of basis states spanning the Hilbert space of quantum many-body 
systems. We test the procedure by working out and analyzing the 
spectral properties of strongly correlated and frustrated quantum spin 
systems. The role and importance of symmetries are investigated. 

PACS numbers: 03.65.-w, 02.70.-c, 68.65.-k, 71.15Nc 

1 Introduction. 

Most microscopic many-body quantum systems are subject to strong inter- 
actions which act between their constituents. The description of some sys- 
tems can be tackled by means of perturbation theory. This may be the case 
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when it is possible to introduce a mean-field concept which is able to include 
quantitatively the main part of the interactions, leaving a remaining residual 
contribution which acts as a more or less small perturbation. 

Often such an approach does not lead to sensible results, in particular 
when one is dealing with realistic quantum spin systems. Non-perturbative 
techniques are needed. During the last decades a considerable amount of 
procedures relying on the renormalization group concept introduced by Wil- 
son [1] have been proposed and tested. Some of them are specifically de- 
vised for quantum spin systems, like the Real Space Renormalization Group 
(RSRG) [3, 4, 6] and the Density Matrix Renormalization Group (DMRG) 
[5, 7, 8]. 

In many cases the study of the spectral properties of quantum systems is 
obtained through the diagonalization of a many-body Hamiltonian in Hilbert 
space spanned by a complete, in general infinite or at least very large set of ba- 
sis states although the information of interest is restricted to the knowledge 
of a few low energy states generally characterized by collective properties. 
Consequently it is necessary to manipulate very large matrices in order to 
extract a reduced quantity of informations. 

Recently we proposed a non-perturbative approach which tackles this 
question [2]. The procedure consists of an algorithm which implements a 
step by step reduction of the size of Hilbert space by means of a projection 
technique. It relies on the renormalization concept following in spirit for- 
mer work based on this concept [9, 10, 11]. Since the reduction procedure 
does not act in ordinary or momentum space but in the zero- dimensional 
Hilbert space like in the procedure developed in ref. [19], it is in principle 
applicable to all types of microscopic quantum systems, in contradistinction 
with methods like the DMRG procedure which works in ordinary space and 
is specifically applicable to Id quantum lattice systems. 

In the present work we implement this algorithm as a preliminary test of 
the practical efficiency and accuracy of the method when applied to strongly 
interacting systems which cannot be treated by means of perturbative meth- 
ods. Second we want to see how far it is able to deliver physical information 
about the properties of the many-body systems it is aimed to describe. We 
use quantum spin systems as first test probes. 
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The outline of the paper is the following. In section 2 we recall the 
essential steps leading to the derivation of the equation which governs the 
evolution of the coupling strength attached to the interaction. The reduction 
formalism is universal in the sense that it works for any kind of many-body 
quantum system. Section 3 is devoted to the application of the algorithm 
to frustrated quantum spin ladders with two legs and one spin per site. We 
analyze the outcome of the applied algorithm on systems of different sizes 
and characterized by different coupling strengths by means of numerical ex- 
amples, with bases of states developed in the SU(2) and S'0(4)-symmetry 
scheme. General conclusions and further planned investigations and devel- 
opments are drawn in section 4. 

2 The reduction algorithm. 

2.1 General concept: the space reduction procedure. 

We consider a system of quantum objects (particles, spins) which are char- 
acterized by a discrete spectrum. The system is governed by a Hamiltonian 



|<?i , g 2 , ' •• , 9p | — ► and acts in a Hilbert space of dimension 

N. The spectrum is obtained from 



<^\^\ N \g( N ^)) , i = 1, • • • , JV| depend on the set of coupling constants 

{g( N ^}- If the relevant quantities of interest are for instance M eigenval- 
ues out of the original set it makes sense to try to define a new effective 
Hamiltonian H^ M \g^ M ^) whose eigenvalues reproduce the M selected states 
and verifies 





(1) 



H^\g^f\g^)) = Hg {M) )\^f\9 {M) )) 



(2) 



with the constraints 



W M) ) = Hg {N) ) 



(3) 
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for i = 1,...,M. If this can be realized Eq. (3) implies a relation between 
the coupling constants in the original and reduced space 

with fc = The effective Hamiltonian H^ M \g^ M ^) may not be rigor- 

ously derivable from H^ N \ It should be constructed so that it optimizes the 
overlap between the original and reduced set of eigenstates. We show next 
how this space reduction may be implemented in practice. 

2.2 Reduction procedure and renormalization of the 
coupling strengths. 

We sketch the procedure which leads from Eq. (1) to Eq. (2). Details can be 
found in ref. [2]. 

We consider a system described by a Hamiltonian depending on a unique 
coupling strength g which can be written as a sum of two terms 

H = H + gH 1 (5) 

The Hilbert space of dimension N is spanned by an a priori arbitrary 
set of basis states i — 1, • • • , N}. An eigenvector 1*^) can be written 

as 

N 

\< ) )-Y.^\9 w )m (6) 

i=i 

where the amplitudes {a[^ \g^)} depend on the value g( N ) of g in «W. 

Using the Feshbach formalism [12] the Hilbert space may be decomposed 
into subspaces by means of the projection operators P and Q, 

H (N) = p H (N) + QH (N) (7) 

In practice the subspace PTiS N ^ is chosen to be of dimension dim P1iS N ^ = 
N — 1 by elimination of one basis state. The projected eigenvector P^^) 
obeys the Schroedinger equation 

H eff (\W)P\*W) = \WP\*W) . (8) 
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where H e ff(\\ ') is the effective Hamiltonian which operates in the sub- 
space PHS N \ It depends on the eigenvalue which is the eigenenergy 
corresponding to l^ 7 ^) in the initial space 1iS N \ The coupling g( N ) which 
characterizes the Hamiltonian in 1iS N ^ is now aimed to be changed into 
such a way that the eigenvalue in the new space TiS N 1 - ) is the same 
as the one in the complete space 

A^ _1) = \[ N) (9) 

The determination of g^ N ~^ by means of the constraint expressed by Eq. (9) 
is the central point of the procedure. In practice the reduction of the vector 
space from iV to N — 1 results in a renormalization of the coupling constant 
from g( N ) to g^ N ~^ preserving the physical eigenenergy A^. 

In the sequel P\ty±*^) is chosen to be the ground state eigenvector and 
= \[ N ^ = Ai the corresponding eigenenergy. In ref. [2] it is shown how 
g^ N ~^ can be obtained as a solution of an algebraic equation of the second 
degree, see Eq. (17). The reduction procedure is iterated in a step by step 
decrease of the dimensions of the vector space, JVKiV-lKi\f-2H... 
leading at each step Ho a coupling strength g( N ~ k ) which can be given as 
the solution of a flow equation in a continuum limit description of the Hilbert 
space [2, 21]. The procedure can be generalized to Hamiltonians depending 
on several coupling constants. 



2.3 Reduction algorithm. 

1— Consider a quantum system described by an Hamiltonian which 
acts in an N- dimensional Hilbert space. 

2— Compute the matrix elements of the Hamiltonian matrix in a def- 
inite basis of states = l,...,iV}. The diagonal matrix elements 
{ej = (&i\H N \$i)} are arranged in increasing order with respect to the {ej}. 

3— Use the Lanczos technique to determine A^ and {^[^(g^)). \[ N ^ may 
be chosen as the experimental value Ai if known. 

4— Fix g^^ 1 ) as described in section 2.2. Take the solution of the algebraic 
second order equation closest to g( N ) (Eq. (17) in [2]). 

5— Construct H^ N ^ = H + g( N ~^ Hi by elimination of the matrix elements 
of H( N ) involving the state |$tv). 
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6— Repeat the procedures 2, 3, 4 and 5 by fixing at each step k X\ ' = 

7— The iterations may be stopped at N = N min corresponding to the limit 
of space dimensions for which the spectrum gets unstable. 

2.4 Preliminary remarks. 

• The procedure is aimed to generate the energies and other physical 
properties of the ground state and low-energy excited states of strongly 
interacting systems. 

• The implementation of the reduction procedure asks for the knowledge 
of Ai and the corresponding eigenvector j 1 !^ - ^) at each step k of the 
reduction process. The eigenvalue Ai is chosen as the physical ground 
state energy of the system. Eigenvalue and eigenvector can be obtained 
by means of the Lanczos algorithm [15, 16, 8] which is particularly well 
adapted to very large vector space dimensions. This algorithm is used 
here in order to fix Ai and |\I f ^ v_ ^). 

• The process does not guarantee a rigorous stability of the eigenvalue Ai. 
Indeed one notices that |\J/ 1 JV_fe which is the eigenvector in the space 
n (N-k-i) and the pro j ec ted state P\^[ N ~ k) ) of \^f{ N ~ k) ) into H^'^ 
may differ from each other. As a consequence it may not be possible 
to keep \[ k ~^ rigorously equal to \[ k ^ = Ai. In practice the degree 
of accuracy depends on the relative size of the eliminated amplitudes 
{ a ik(9^)}- This point will be tested by means of numerical estima- 
tions and further discussed below. 

• The algorithm and different points which have been quoted above will 
now be developed in applications of the procedure to explicit models, 
here frustrated spin ladders. The applicability of the algorithm will be 
tested in different symmetry schemes. 

• The Hamiltonians of the considered ladder systems are characterized by 
a fixed total magnetic magnetization M tot . We shall work in subspaces 
which correspond to fixed M tot . The total spin S to t is also a good 
quantum number which defines smaller subspaces for fixed M tot . We 
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li 2i 3i 4i 5i ii s Jl= l/2 L 




1 2 ^ 2 2 3 2 4 2 5 2 ; 2 s h L 



1 2 3 4 5 / L 

Figure 1: Top: the original spin ladder. The coupling strengths are indicated 
as given in the text. Bottom: The ladder in the SO(4) representation. See 
the text. 

shall not introduce them here because projection procedures on Stot are 
time consuming. Furthermore we want to test the algorithm in large 
enough spaces although not necessarily the largest possible ones in this 
preliminary tests considered here. 

3 Application to frustrated two-leg quantum 
spin ladders. 

3.1 The model 

3.1.1 SU(2)-symmetry framework. 

Consider spin-1/2 ladders [13, 14] described by Hamiltonians of the following 
type and shown in Fig. 1. 

L 

H (s,s) = Jt y^ Sii g i2 + Jj s^Sj! + Ji^2 Sj 2 s h + Ji c ^2s h s h (10) 

i=l <ij> <ij> (ij) 

(ij) 

The indices 1 or 2 label the spin 1/2 vector operators s ik acting on the 
sites i on both ends of a rung, in the second and third term i and j label 
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nearest neighbours, here j = i + 1 along the legs of the ladder. The fourth 
and fifth term correspond to diagonal interactions between sites located on 
different legs, j — % + 1. L is the number of sites on a leg (Fig. 1) where 
Jic — -he = J c - The coupling strengths J t , Ji, J c are positive. 

As stated above the renormalization is restricted to a unique coupling 
strength, see Eq. (5). It is implemented here by putting H = and = 
g( N ^H 1 where = J t and 

L 

Hi = ^2 Si i Si 2 + 7ti Yl ( Si i^'i + + 7c Yl ( Si i s i* + ^SjJ • (11) 

i=l <ij> <ij> 

where 7^ = J;/Jj, 7 C = J c / Jt- These quantities are kept constant and 
g( N ) — j t w jii De subject to renormalizationin the reduction process. 

The basis of states {\&k) k — 1, . . . , N} is chosen as 

2L 

|$ fc ) = |l/2 mi, ...,1/2 m^, ...,1/2 m 2L , ^ m, = M tot = 0) 

i=i 

with {mi = +1/2,-1/2}. 

3.1.2 SO(4)-symmetry framework. 

The basis of states may be written in an S'0(4)-symmetry scheme. Differ- 
ent choices of bases may induce a more or less efficient reduction procedure 
depending on the strength of the coupling constants J t , Ji, J c . This point is 
investigated here. 

By means of a spin rotation [17, 18] 

s ll = ^(S l + R t ) . (12) 

s i2 = \{Si - Ri) . (13) 
the Hamiltonian Eq.(10) can be expressed in the form 



8 



Hi s, R) = J± J2(S> - R>) + Jl J2 S t S, + J 2 J2 RiRj (14) 

i=l <ij> <ij> 

The structure of the corresponding system is shown in the lower part of Fig. 
I. Here J x = (J t + J c )/2, J 2 = (Jj — J c )/2 and as before J lc = J 2c = J c . The 
components S^ + \ S^~\ and R\ , R^ of the vector operators Si and 
Ri are the 5*0(4) group generators and < ij > denotes nearest neighbour 
indices 



KH)(10) . v (10)(1-1)n 



(-)* 



(z) = X (U)(U) _ x (l-l)(l-l) 



R^ = V2{x\ ll)m - X 



(11)(00) „(00)(1-1)n _ „(-) 



)=R\ 



where 



to _^(io)(oo) +x (oo)(io)^ 



In this framework the states {\SiMi)} are defined as 

\SiMi) = (V2 mi 1/2 mal^M,)!!^ mi)i|l/2 m 2 ) 2 



mi,m2 



along a rung are coupled to Si — or Si — 1. Spectra are constructed in this 
representation as well as in the SU (2) representation and the basis of states 
{|$fc)} takes the form 



|$ fc ) = I^Mx, S,M t , 5 L M L , ^ Mj = M tot = 0) 



i=l 
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3.2 Test observables 

In order to quantify the accuracy of the procedure we introduce different 
test quantities in order to estimate quantitatively deviations between ground 
state and low excited state energies in Hilbert spaces of different dimensions. 
The stability of low-lying states can be estimated by means of 

/ (AO _ (AT-fch 

p(i) = |lJ —i q x 100 with i=l,..., 4 (15) 

4 

where e\ N ~ k ^ = \\ N ~ k ^/2L corresponds to the energy per site at the ith phys- 
ical state starting from the ground state at the fcth iteration in Hilbert space. 
These quantities provide a percentage of loss of accuracy of the eigenenergies 
in the different reduced spaces. 

A global characterization of the ground state wavefunction in different 
symmetry schemes can also be given by the entropy per site in a space of 
dimension n 

* = -^£Wi with P* = \& n) \^r ) )\ 2 = \aV\ 2 (16) 
i=i 

which works as a global measure of the distribution of the amplitudes {a^} 
in the physical ground state. 



3.3 Spectra in the SU(2)-symmetry framework. 

We apply the reduction algorithm to ladders with two legs, different num- 
bers of sites and different values of the coupling strengths. Results obtained 
with an SU (2)-symmetry basis of states are shown in Figs. (2 - 4) and Fig. 6. 

3.3.1 First case: L= 6, J t =15, J t =5, J c =3 

We choose the basis states in the framework of the M-scheme corresponding 
to subspaces with fixed values of the projection of the spin of the 
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M tot = 0. 



In the present case J t > Ji, J c . The dimension of the subspace is reduced 
step by step as explained above starting from N = 924. As stated in section 
2.3 the basis states are ordered with increasing energy of their diago- 

nal matrix elements {e^} and eliminated starting from the state with largest 
energy, e N . 

As seen in Fig. (2a) the ground state of the system stays stable down to 
n ~ 50 where n is the dimension of the reduced space. The coupling constant 
J t does not move either up to n ~ 300. Figs.(2a-b) show the evolution of the 
first excited states which follows the same trend as the ground state. Devia- 
tions from their initial value at N = 924 can be seen in Fig.(2c-d) where the 
p(i)'s defined above represent these deviations in terms of percentages. 

For n < 50 the spectrum gets unstable, the renormalization of the cou- 
pling constant can no longer correct for the energy of the lowest state. Indeed 
the coupling constant Jt increases drastically as seen in Fig.(2e). The reason 
for this behaviour can be found in the fact that at this stage the algorithm 
eliminates states which have an essential component in the state of lowest 
energy. The same message can be read on Fig.(2f), the drop in the entropy 
per site s is due to the elimination of sizable amplitudes {an}. 

3.3.2 Second case: L= 6, J t =5.5, J;=5, J c =3 

Contrary to the former case the coupling constant J t along rungs is now of 
the same strength as Ji, J c . Results are shown in Fig. (3). The lowest energy 
state is now stable down to n ~ 100. This is also reflected in the behaviour 
of the excited states which move appreciably for n < 200. Fig.(3e) shows 
that the coupling constant J t starts to increase sharply between n = 300 and 
n = 200. It is able to stabilize the excited states up to about n = 200 and 
the ground state up to n = 70. The instability for n < 70 reflects in the 
evolution of the p(i)'s, Figs.(3c-d) which get of the order of a few percents. 
The entropy Fig.(3f) follows the same trend. 

Comparing the two cases above and particularly the entropies Fig.(3f) 
and Fig.(2f) one sees that the stronger J t the more the amplitude strength of 
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the ground state wavefunction is concentrated in a smaller number of basis 
state components. The elimination of sizable components of the wavefunction 
leads to deviations which can be controlled up to a certain limit by means 
of the renormalization of J t . One sees that large values of J t favour a low 
number of significative components in the low energy part of the spectrum 
in a £77(2) symmetry framework. 

A confirmation of this trend can be observed in Figs.(4a-f) where J t = 2.5. 
The rates of destabilisation of the excited states are higher than in the for- 
mer cases as it can be seen in Figs.(4c-d). This point is also reflected in the 
behaviour of the entropy s which is larger than in the former case for n = N 
and decreases more rapidly with decreasing n, Fig.(4f). 

3.3.3 Third case: L= 8, J t =15, J/=5, J c =3 

For L = 8 the Hilbert space is spanned by iV = 12870 basis states with 
M to t = 0. The results are shown in Figs.(6a-f). The stability of the spectrum 
with decreasing space dimension is stronger than the stability observed for 
7 = 6. Indeed iin/N defines the ratio of the number of states in the reduced 
space over the total number of states one finds p(l) ~ 0.8% and p(2) ~ 0.8% 
when n/N ~ 0.07 for 7 = 6. For L = 8 p(l) ~ 0.8% when n/N ~ 0.007 and 
p(2) ~ 0.5% when n/N ~ 0.02. This shows a sizable improvement in the 
stability of the spectrum, at least in the specific domain where the coupling 
strength J t is large compared to the others. 

The evolution of the spectrum and its stability with decreasing J t follows 
the same trend as in the case where L = 6. 



3.3.4 Remarks 

In Fig. (2a) it is seen that the ground state shows "bunches" of energy fluc- 
tuations. The peaks are intermittent, they appear and disappear during the 
space dimension reduction process. They are small in the case where J t = 15 
but can grow with decreasing J t as it can be observed for J t = 2.5. The 
subsequent stabilization of the ground state energy following such a bunch 
shows the effectiveness of the coupling constant renormalization which acts 
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in a progressively reduced and hence incomplete basis of states. 

These bunches of fluctuations are correlated with the change of the num- 
ber of relevant amplitudes (i.e. amplitudes larger than some value e as ex- 
plained in the caption of Fig. (5)) during the reduction process. 

Consider first the case where J t > Ji, J c - One notices in the caption of 
Fig. (5a) that up to n ~ 300 the number of relevant amplitudes defined in 
Fig. (5) stays stable like the ratios {p(i)} in Figs.(2c-d). For 158 < n < 300 
these ratios change quickly. A bunch of fluctuations appears in this domain 
of values of n as seen in Figs.(2c-d) and correspondingly the number of rel- 
evant amplitudes decreases steeply. For 60 < n < 158 the ratios {p(i)} stay 
again stable as well as the number of relevant amplitudes. The {p{i)} in 
Fig.(2c-d) almost decrease back to their initial values. The same explanation 
is valid for L = 8 in Fig. (6). The analysis shows that these bunches of fluctu- 
ations signal the local elimination of relevant contributions of basis states to 
the physical states in the spectrum. The stabilization of the spectra which 
follows during the elimination process shows that renormalization is able to 
cure these effects. 

In the case where J t < Ji, J c shown in Fig. (5b) the relevant and irrel- 
evant amplitudes move continuously during the reduction process and the 
corresponding {p(i)} do no longer decrease to the values they showed before 
the appearance of the bunch of energy fluctuations as seen in Figs.(6c-d). It 
signals the fact that the coupling renormalization is no longer able to com- 
pensate for the reduction of the Hilbert space dimensions. 



3.4 Spectra in the SO(4)-symmetry framework 

The reduction algorithm is now applied to the system described by the Hamil- 
tonian H^ S ' R ^ given by Eq. (14) with a basis of states written in the SO (A) 
symmetry framework. Like above we consider two cases corresponding to 
large and small values of J t relative to the strengths of the other coupling 
parameters. 
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3.4.1 Reduction test for L = 6, J t = 15, 2.5 and J;=5, J c =3 

Figs. (7) show the behaviour of the spectrum for a system of size L = 6. A 
large value of J t , (J t = 15), favours the dimer structure along rungs in the 
lowest energy state and stabilizes the spectrum down to small Hilbert space 
dimensions. This effect is clearly seen in Fig. (7a), the ground state is very 
stable. The excited states are more affected, see Figs. (7b), although they 
do not move significantly, Figs.(7c-d). The renormalization of the coupling 
strength J t starts to work for n ~ 50. 

The situation changes progressively with decreasing values of J t . Figs. (8) 
show the case where Jt = 2.5. The ground state energy experiences sizable 
bunches of fluctuations like in the 577(2) scheme, but much stronger than 
in this last case. The same is true for the excited states which is reflected 
through all the quantities shown in Figs. (8), in particular J t) Fig.(8e). The 
arguments used in the 5£/(2)-scheme about relevant and irrelevant ampli- 
tudes are also valid here. 

The result shows that the renormalization procedure is quite sensitive to 
the symmetry scheme chosen in Hilbert space. It is expected that essential 
components of the ground state wavefunction get eliminated early during 
the process when the rung coupling gets of the order of magnitude or smaller 
than the other coupling strengths. 



3.5 Summary 

The present results lead to two correlated remarks. The efficiency of the 
algorithm is different in different sectors of the coupling parameter space. 
In the case of the frustrated ladders considered here the algorithm is the 
more efficient the stronger the coupling between rung sites J t . Second, this 
behaviour is strongly related to the symmetry representation in which the 
basis of states is defined. The 577(2) representation leads to a structure of 
the wavefunctions (i.e. the size of the amplitudes of the basis states) which is 
very different from the one obtained in the 5*0(4) representation. For large 
values of J t the spectrum is more stable in the SO (4) scheme. For small 
values of Jt the stability is better realized in the 577(2) scheme. Finally, in 
the regime where J t > Ji, J c , one observes that the reduction procedure is the 
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more efficient the closer Ji to J c . This effect can be understood and related 
to previous analytical work in the 5*0(4) framework [20]. 



4 Conclusions and outlook. 

In the present work we tested and analysed the outcome of an algorithm 
which aims to reduce the dimensions of the Hilbert space of states describing 
strongly interacting systems. The reduction is compensated by the renor- 
malization of the coupling strengths which enter the Hamiltonians of the 
systems. By construction the algorithm works in any space dimension and 
may be applied to the study of any microscopic iV-body quantum system. 
The robustness of the algorithm has been applied to frustrated quantum spin 
ladders. 

The analysis of the numerical results obtained in applications to quantum 
spin ladders leads to the following conclusions. 

• The stability of the low-lying states of the spectrum in the course of 
the reduction procedure depends on the relative values of the coupling 
strengths. The ladder favours a dimer structure along the rungs, i.e. 
stability is the better the larger the transverse coupling strength J t . 

• The evolution of the spectrum depends on the initial size of Hilbert 
space. The larger the initial space the larger the ratio between the 
initial number of states and the number of states corresponding to the 
limit of stability of the spectrum. 

• The efficiency of the reduction procedure depends on the symmetry 
frame in which the basis of states is defined. It appears clearly that the 
evolution of the spectrum described in an SU(2) scheme is significantly 
different from the evolution in an SO (A) scheme. This is again under- 
standable since different symmetry schemes partition Hilbert space in 
different ways and favour one or the other symmetry depending on the 
relative strengths of the coupling constants. 
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Local spectral instabilities appearing in the course of the reduction pro- 
cedure are correlated with the elimination of basis states with sizable 
amplitudes in the ground state wavefunction. One or another repre- 
sentation can be more efficient for a given set of coupling parameters 
because it leads to physical states in which the weight on the basis 
states is concentrated in a different number of components. This point 
is strongly related to the correlation between quantum entanglement 
and symmetry properties which are presently under intensive scrutiny, 
see f.i. [22] and refs. therein. 

Further points are worthwhile to be investigated: 

• In the present approach the sequential reduction of space dimensions 
followed an energy criterion. It might be judicious to classify the se- 
quence of states to be eliminated starting with those which have the 
smallest amplitude in the ground state wavefunction. The two proce- 
dures should however be correlated if not equivalent. 

• We expect to extend the study to systems of higher space dimensions, 
f.i. 2d. 

• The present approach relies on an algorithm which is able to recognize 
the existence of first and higher order critical points [2]. It is of in- 
terest to apply the algorithm in the neighbourhood of such points. Its 
behaviour could help to identify them. Work on this point is under way. 

• The algorithm can be extended to systems at finite temperature [21] 
and more than one coupling constant renormalization. 

The authors would like to thank Dr. A. Honecker for interesting discus- 
sions. One of us (T.K.) acknowledges the help of Drs. E. Caurier and F. 
Nowacki on technical aspects concerning the implementation of the Lanczos 
algorithm. 
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Figure 2: SU(2) — scheme. N is Hilbert space dimension. The {ei,i 
1,2,3,4} are the energies of the ground and excited states per site. L = 
sites along a leg. J t = 15, Ji = 5, J c = 3 
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Figure 3: SU(2) — scheme. N is Hilbert space dimension. The {e^, i 
1,2,3,4} are the energies of the ground and excited states per site. L = 
sites along a leg. J t = 5.5, J[ = 5, J c = 3 
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Figure 4: SU(2) — scheme. N is Hilbert space dimension. The {e^, i = 
1, 2, 3, 4} are the energies per site. L = 6 sites along a leg. J t = 2.5, Ji = 5, 

Jr =3 
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Figure 6: SU(2) — scheme. N is Hilbert space dimension. The {e^ 
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Figure 7: SO(4) — scheme. N is Hilbert space dimension. The {ei,i = 
1,2,3,4} are the energies per site. L = 6 sites along the chain. J t = 15, 
Jl = 5, J r = 3 
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Figure 8: 5*0(4) — scheme. N is Hilbert space dimension. The {ei,i = 
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